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In the present paper we study the behavior of the normalized I-Q relation for neutron stars in a 
particular class of f(R) theories of gravity, namely the R 2 gravity that is one of the most natural and 
simplest extensions of general relativity in the strong field regime. We study both the slowly and 
rapidly rotating cases. The results show that the I-Q relation remain nearly equation of state inde¬ 
pendent for fixed values of the normalized rotational parameter, but the deviations from universality 
can be a little bit larger compared to the general relativistic case. What is the most interesting in our 
studies, is that the differences with the pure Einstein's theory can be large reaching above 20%. This 
is qualitative different from the majority of alternative theories of gravity, where the normalized I-Q 
relations are almost indistinguishable from the general relativistic case, and can lead to observational 
constraints on the f(R) theories in the future. 

PACS numbers: 


I. INTRODUCTION 

Studies of generalized theories of gravity are becoming more and more intense in the last decade. There are both 
theoretical and observational motivations for this. On one hand the theories trying to unify all the interactions 
predict that the standard Einstein-Hilbert action should be modified. On the other hand it was shown in many 
cases that studying generalizations of Einstein's gravity can give us a deeper understanding of general relativity 
(GR) itself. On the observational front still remain phenomena that do not fit very well in the standard framework, 
such as the accelerated expansion of the universe, and that is why modifications of the theory of gravity are often 
employed as an alternative explanations. One should keep also in mind, that even though general relativity is very 
well tested in the weak field regime, the strong field remains essentially unconstrained that leaves space for a variety 
of modifications. 

One of the most natural generalizations of Einstein's theory of gravity are the f(R) theories, where the Ricci scalar 
R in the Einstein-Hilbert action is replaced by some function of R. Such modification has a theoretical motivation for 
example from the quantum field theory in curved spacetime. f(R) theories are also widely used as an alternative 
explanation of the dark energy phenomena which places them amongst the most popular and widely explored 
alternative theories of gravity. Most of the studies on f(R) theories though are in cosmological aspect in relation 
to the accelerated expansion of the universe fl}[3|. The examinations of the astrophysical manifestations of these 
theories is more scarce and this would be the main focus of our paper. 

Natural objects to study within the f(R) theories of gravity at astrophysical scales are the neutron stars, where 
the strong gravity effects are non-negligible. The neutron stars within the f(R) theories can differ significantly from 
their GR counterpart [5H6j which makes them a very good candidate to test f(R) theories on astrophysical scales. 
Unfortunately one has to pay a high price - the nuclear matter equation of state (EOS) at densities as high as the ones 
in the neutron star cores, is still unknown. That is why in many cases the deviations coming from the generalizations 
of Einstein's gravity are comparable or even smaller than the deviations resulting from the uncertainties in the EOS. 

A way to circumvent this problem is to search for predictions or derive relations that are independent of the EOS. 
This was exactly the idea in 13(81 where the famous f-Love-Q relations were discovered. These relations connect the 
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normalized neutron star moment of inertia I, quadrupole moment Q and the tidal Love number A. They were later 
explored in the regime of rapid rotation l!9tiTTll . large tidal deformations Hl2 l, in the presence of magnetic fields Hi 31 . 
for higher multipole moments [14] and in the presence of anisotropic pressure [15]. We should note that also a variety 
of other (nearly) EOS independent relations exists in the literature connecting different neutron star properties, their 
oscillations spectrum, etc. (see for example 1IT6H26I ). 

The most important property of the f-Love-Q relations is that they are practically independent of the EOS (for 
moderate magnetic fields). Several applications were proposed and one of the most important is braking the de¬ 
generacy between the spins and the quadrupole moment of neutron star inspirals, and testing alternative theories 
of gravity. In connection to the latter, the I-Love-Q relations were examined in a variety of alternative theories of 
gravity, such as the dynamical Chern-Simons gravity f7J (SJ, Eddington-inspired Born-Infeld gravity l27l , Einstein- 
Gauss-Bonnet-dilaton theory l28l and scalar-tensor theories of gravity l 29]|32| . In all of these theories the f-Love-Q 
relations are pretty much EOS independent. With the exception of the dynamical Chern-Simons gravity, the resulting 
relations are quite similar to the GR case and they can not be used to test the alternative theories of gravity 1 . 

The purpose of the present paper is to explore the normalized I-Q relations for f(R) theories of gravity and more 
specifically for the R 2 gravity, where f(R) — R + aR 2 . We will examine both the slowly and rapidly rotating regime. 
The main goals is to determine if the relations are still EOS independent and whether significant deviations from GR 
can be observed. We will use the fact that /(R) theories are mathematically equivalent to a particular class of scalar- 
tensor theories with a nonzero potential of the scalar field. As we will demonstrate below, exactly the presence of 
such nonzero potential leads to interesting results and makes the problem (and the conclusions) qualitative different 
from the one considered in 1501 . 


II. BASIC EQUATIONS 


In our calculations we will use the fact that f(R) theories are mathematically equivalent to a specific class of scalar- 
tensor theories. Here we will give very briefly the main points and we refer the reader to [4]46| where the problem is 
discussed in detail. 

The f(R) gravity action can be written in the following general form 

s — J d X \J—Sf(R) T S ma tter(gf(V/ X)> (1) 

where R is the Ricci scalar curvature with respect to the space-time metric g^ v , S ma tter is the action of the matter, and 
the matter fields are denoted by This action is mathematically equivalent to the following scalar-tensor theory 
(STT) action 

S = J [R* - 2 g*l lv d }l <pd v tp -V(cp)] + S matter (A 2 {cp)g* v ,x), (2) 

where the coupling function A(cp') and the scalar-field potential V (<p) have the following form 

A 2 (?) - HK V(cp) - A\<p) (r^ - /(R)) . (3) 


Here R* is the Ricci scalar curvature with respect to the metric g* lv and cp is the scalar field. In this paper we will 
concentrate on the so-called R 2 gravity defined by /(R) = R + aR 2 . In this case the scalar-field potential takes the 
form 


1 / _ 2£\ 2 


(4) 


1 One should keep in mind that this is true only for the normalized relations. The unnormalized quantities can deviate significantly and poten¬ 
tially lead to some observational effects. 
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The STT action § is written in the so-called Einstein frame, that is not the physical one and it is introduced because 
it simplifies the field equations substantially. The physical quantities, such as mass, distance, etc., are measured in 
the Jordan frame, where the two frame are connected by a conformal transformation of the metric and a redefinition 
of the scalar field. We will not go into details and we refer the reader to J4| |331 for an extensive discussion of the 
problem. 

We will consider stationary, axisymmetric and asymptotically flat solutions of the field equations describing rotat¬ 
ing compact starts. Thus the Einstein frame spacetime metric can be presented in the following form 

dsl — —e 2v dt 2 + p 2 B 2 e~ 2v (dcp — codt) 2 + e 2 ^~ 2v (dp 2 + dz 2 ), (5) 

where all the metric functions depend on the coordinates p and z only. Let us comment on the asymptotic behavior 
of the metric functions and the scalar field. It is derived from the reduced field equations that can be found in 0. It 
is more convenient to present the asymptotic behaviour in the quasi-isotropic coordinates r and 0 defined by 


p — r sin 6, z — r cos 9. 

Using these coordinates and keeping only terms up to the order of r~ 3 , one can obtain 
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where M and / are the mass and the angular momentum, b, v 2 an d <p 2 are constants and P 2 (cos0) is the second 
Legendre polynomial. The scalar field on the other hand decreases exponentially at infinity and one can show that 


V 





( 11 ) 


for large values of r, where m ( p is the scalar field mass equal to m (/ , — 1/ (\/6a). One can make two important 
conclusions from eq. ©• First, the scalar charge of the scalar field is zero 2 and therefore it does not contribute to the 
asymptotic of the metric functions (7| - (10) contrary to the case of the scalar-tensor theories (with massless scalar 
field) admitting scalarization f30]]. Second, the exponential decay of the scalar field is controlled by the parameter a 
and smaller a corresponds to more rapid decay. Thus the general relativistic case corresponds to the limit a —> 0. 

Since the scalar field does not contribute to the asymptotic of the metric, the quadrupole moment would have the 
same form as in general relativity. Thus one can obtain Il34ll35l 


Q = -v 2 - 


4 

3 



( 12 ) 


In our calculations we will use also the moment of inertia I defined in the usual way 



with O being the angular velocity of the star. 


( 13 ) 


2 


The scalar charge in usually defined as the coefficient in front of the 1/r term in the scalar field expansion at infinity. 
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Equations (12) and (13) are written in the Einstein frame. Therefore, one has to transform them in the physical 
Jordan frame. But, taking into account the exponential decay of scalar field, one can easily show that the moment of 
inertia and the quadrupole moment are the same in both frames for the particular scalar-tensor theory and therefore 
for the R 2 gravity we consider in the present paper. Therefore, formulae and |l3| are valid also in the physical 
Jordan frame. 

In our calculations we use the dimensionless parameter a —> a/R jy where Rq is one half of the solar gravitational 
radius Rq — 1.47664 km (i.e. the solar mass in geometrical units). 


III. NUMERICAL RESULTS 

In our calculations we employed four equations of state (EOS) that cover a wide range of stiffness - APR, SLy4, 
FPS and Shen2D (the zero temperature limit of the Shen EOS). The first two EOS are considered as standard since 
they are above the two solar mass barrier and have radii in the preferred range according to the observations 13614391 . 
The FPS EOS is softer with maximum mass below the 2 Mq . The Shen2D EOS on the other hand is stiff and has high 
maximum mass, but the typical radii are around 14-15km that is larger compared to the observational estimates. 
Even though EOS FPS and Shen2D do not fit all the observations well, we have chosen them in order to cover a 
larger range of stiffness and to check up to what extend our relations are EOS independent. The rotating neutron 
star solutions in R 2 gravity are obtained using an extended version of the rns code 0. 

The goal of the present paper is to study the relation between the moment of inertia I and the quadrupole moment 
Q of neutron stars in f(R) theory in connection to the famous I-Love-Q relations. In order to obtain (nearly) equation 
of state independent relations, I and Q have to be properly normalized. A natural choice is I = 1/ M 3 and Q = 
Q/ (M 3 y 2 ), where X = I / Ad 2 . Since we are considering rotating models, we will have two main goals - first is to 
determine what the deviations from pure general relativity are, and second - to check up to what extend the relations 
are equation of state independent. We should note that a different normalization might make the relations more or 
less universal and increase or decrease the differences between the gravitational theories. A conjecture was made in 
[40J that always exists a normalization that can make a given relation universal. But one should not choose aways the 
most universal normalization, but instead the most relevant one from both physical and observational point of view. 
The particular normalization in the current paper is the standard one used in all the previous papers on I-Love-Q 
relations, including relations in alternative theories of gravity. It has a number of physical justifications and makes 
the comparison of the results much easier. 

We work with two different rotational parameters that are popular in the literature - the normalized rotational 
frequency vM (v is the rotational frequency of a star in kHz, M is the mass in solar masses) and the normalized 
angular momentum / / M 2 (// M 2 is dimensionless in units c — G — 1). The I — Q relations for sequences of neutron 
star models with constant value of vM are presented in the left panel of Fig. [ljand with constant //M 2 - in the right 
panel of Fig. [T] The sequences range from the slow rotation limit (vM — 0.3 and //M 2 = 0.1) to very fast rotation 
(vM — 2.7 and //M 2 = 0.6). For every single vM three cases are shown - the general relativistic limit (a — 0) and R 2 
gravity with a = 2 and a = 10 4 (a = 10 4 gives nearly the maximum possible deviation from GR). In the case of //M 2 
we show the results only for a — 0 and a — 10 4 in order to have a better visibility. 

In the bottom part of the plots in Fig. [Tjthe deviations A I of the /(R) gravity results from the GR case are shown, 
defined as 

AI = |rGR -~ lFR| . (14) 

Igr 

Here Igr and JpR are the corresponding values of the moment of inertial in the GR and /(R) cases respectively and 
the plotted values of AI are in percentages. 

The most important result is that AI reaches above 20% for the less massive models (with larger Q)/ slower ro¬ 
tation and large values of a. Such deviation is much higher that the one observed for other alternative theories of 
gravity such as scalar-tensor theories which admit scalarization l30| . This is surprising at first sight because, as we 
described above, we are using the mathematical equivalence between /(R) theories and scalar-tensor theories in 
our calculations. But there is a very important difference between the particular class of STT considered here that 
is equivalent to the /(R) theories (see eqs. and the STT considered in |[30jl33| . It steams from the fact that the 
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FIG. 1: I — Q relation for sequences with fixed rotational parameter vM (right panel) and fixed //M 2 (left panel). The lower panel 
for both figures represents the relative deviation of the results in f(R) gravity from the GR case. 


scalar-field mass in our case is nonzero which leads to its exponential decay (see eq. Thus the scalar field has 

a zero scalar charge and it does not contribute directly to the asymptotic of the metric functions at infinity and the 
quadrupole moment. Therefore eqs. (7| and l |12| are valid both in GR and R 2 gravity. 

This is not the case with the STT considered in Ii30ll33) , which admits spontaneous scalarization and has a massless 
scalar field. A nonzero scalar charge is present there which enters explicitly in the equation for the quadrupole 
moment and roughly speaking it compensates the deviations from GR in the I — Q relation 3 . Such "compensation" 
of the scalar field effect does not exist in the STT considered in this paper (and therefore in f(R) gravity) and that is 
why large differences from GR can be observed. 

The deviations from EOS universality are small, below roughly 1%, for the general relativistic case (a=0) and 
for the case of very large a. The deviations can increase for small a reaching roughly 3%, but this case is not so 
interesting since it gives small differences with GR. Another point is that deviations as large as 3% are observed only 
for more extreme EOS, such as the Shen2D EOS. For the "standard" modern realistic EOS the deviations from EOS 
universality are still within roughly 1%. 

Let us comment on the effect of rotation. From the bottom panel in Fig. [I] it is evident that for fixed values of Q 
the deviations from GR are comparable for different values of the normalized rotational parameters vM and // M 2 . 


3 We should note that only the normalized I — Q relation is close to GR for the STT considered in 1301 . The unnormalized quantities I and Q can 
reach very large deviations from GR especially in the rapidly rotating case 
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Therefore, large differences can be observed already for slow rotation that is very important since a big portion of 
the observed systems have rotational frequencies up to a few hundred Hz. If the sequences with fixed normalized 
rotational parameter terminate at smaller values of Q, i.e. at large masses, the deviations from GR are relatively small 
that naturally come from the fact that they are close to the black hole limit where the I — Q relation does not depend 
on the internal structure of the star (this is observed in Fig. [ljespecially in the case of fixed vM). The deviations from 
EOS universality are not sensitive to the particular value of vM or //M 2 , i.e. they are of the same order for all of the 
sequences. 

Another important point one can notice is that the effect of rotation is qualitatively the same as the effect of 
increasing the parameter a. Therefore, in order to test the deviations from general relativity we should either have 
additional information about the rotational rate of the star from independent observations, or deal with objects that 
are in general slowly rotating (the I — Q dependences are almost indistinguishable for frequencies less than a few 
hundred Hz). The opposite line of conclusions is also possible - it would be very difficult to test any deviations from 
slow rotation via the f-Love-Q relations since certain modifications of gravity lead to similar results. 


IV. CONCLUSIONS 

In the present paper we have studied the normalized I — Q dependence for rotating neutron stars in R 2 theory of 
gravity. Both the slowly and the rapidly rotating cases are examined. In our calculation we made use of the fact that 
the f(R) theories are mathematically equivalent to a particular class of scalar-tensor theories with a nonzero potential 
of the scalar field. The effective mass of the scalar field in this scalar-tensor representation is backproportional to the 
R 2 gravity parameter a and the scalar field decays exponentially at infinity. Therefore the scalar charge is zero and 
the asymptotic of the metric at infinity and the quadrupole moment have the same form as in pure general relativity 
(contrary to other classes of STT with massless scalar field 1301 1. The rotating neutron star solutions were obtained 
with an extended version of the rns code |5|. 

The most important conclusion from our results is that the dependences can deviate significantly from the GR case 
and the differences can reach above 20% that can be observationally relevant. The deviations are in general larger 
for larger values of the R 2 gravity parameter a, smaller neutron star masses and slower rotation. This is qualitatively 
different from the case of scalar-tensor theories with vanishing potential, where the normalized I — Q relation differs 
only marginally from GR even in the rapidly rotating case. A rough explanation of this phenomena comes from 
the fact that the scalar charge is zero in the scalar-tensor representation of the R 2 gravity and thus it has no explicit 
contribution to the quadrupole formula. The effect of R 2 gravity is qualitatively similar to the effect of changing the 
normalized rotational parameter vM or //M 2 . So one should be careful when trying to distinguish between the two 
effects in the future observations. 

The I — Q relations in R 2 gravity remain nearly EOS independent for fixed values of the normalized rotational 
parameters vM and // M 2 similar to the pure Einstein's theory. The deviations are typically of the order of 1% for 
large values of a or large masses and they increase up to roughly 3% with the decrease of a and M. This is certainly 
below the expected observational accuracy. Also deviations as large as 3% are observed only for more extreme EOS, 
and therefore they are of limited interest. 

It is interesting to note that the considered /(R) theories are one of the few generalized theories of gravity that 
lead to large deviations from GR as far as the I — Q relations are concerned. As we commented above, this is due to 
the fact that the scalar field in the scalar-tensor representation of the R 2 gravity is massive (i.e. has a finite range). 
Hence we expect that other classes of massive scalar-tensor theories will have similar behavior and such a study is 
underway. 

At the end we will briefly comment on the observational perspectives. As we have already pointed out, larger 
deviations in the I — Q relations are observed for lower masses and slow rotation. But this is exactly the case of 
merging neutron stars that are supposed to have typical masses around 1.4Mq and to be slowly rotating due to 
the long evolution time. Probably the /(R) gravity will also change the gravitational wave signal emitted by such 
inspiraling neutron stars similar to the case of scalar-tensor theories l[4lij44ll but additional studies in this direction 
are needed. One of the most promising application of the f-Love-Q relations is in the parameter extraction of such 
merging binaries and therefore one should keep in mind the possible effects that might come from the alternative 
theories of gravity and the /(R) theories in particular. 
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As far as isolated and binary NS are concerned, the precise observation of two of the quantities in the t-Love-Q 
relation would serve as a direct test of GR and can help us to impose constraints on f(R) theories. But we should 
have at least some information about the rotational frequency of these stars and whether they are slowly or rapidly 
rotating, since the deviations coming from Rr gravity are qualitative the same as the effect of rotation. 


Acknowledgements 

DD would like to thank the European Social Fund and the Ministry Of Science, Research and the Arts Baden- 
Wurttemberg for the support. SY and KS would like to thank the Research Group Linkage Programme of the 
Alexander von Humboldt Foundation for the support. The support by the Bulgarian NSF Grant DFNI T02 / 6, Sofia 
University Research Fund under Grant 70 /2015 and "New-CompStar" COST Action MP1304 is gratefully acknowl¬ 
edged. 


[1] T. P. Sotiriou and V. Faraoni, Rev.Mod.Phys. 82, 451 (2010). 

[2] S. Nojiri and S. D. Odintsov, Phys.Rept. 505, 59 (2011). 

[3] S. Capozziello and M. de Laurentis, Phys. Rep. 509,167 (2011), 1108.6266. 

[4] S. S. Yazadjiev, D. D. Doneva, K. D. Kokkotas, and K. V. Staykov, JCAP 1406, 003 (2014). 

[5] S. S. Yazadjiev, D. D. Doneva, and K. D. Kokkotas, Phys. Rev. D 91, 084018 (2015), 1501.04591. 

[6] K. V. Staykov, D. D. Doneva, S. S. Yazadjiev, and K. D. Kokkotas, /. Cosmology Astropart. Phys. 10, 006 (2014), 1407.2180. 

[7] K. Yagi and N. Yunes, Science 341, 365 (2013). 

[8] K. Yagi and N. Yunes, Phys. Rev. D 88, 023009 (2013). 

[9] D. D. Doneva, S. S. Yazadjiev, N. Stergioulas, and K. D. Kokkotas, Ap] 781, L6 (2014), 1310.7436. 

[10] G. Pappas and T. A. Apostolatos, Physical Review Letters 112,121101 (2014), 1311.5508. 

[11] S. Chakrabarti, T. Delsate, N. Giirlebeck, and J. Steinhoff, Physical Review Letters 112, 201102 (2014), 1311.6509. 

[12] A. Maselli, V. Cardoso, V. Ferrari, L. Gualtieri, and P. Pani, Phys. Rev. D 88, 023007 (2013). 

[13] B. Haskell, R. Ciolfi, F. Pannarale, and L. Rezzolla, ArXiv e-prints (2013). 

[14] K. Yagi, K. Kyutoku, G. Pappas, N. Yunes, and T. A. Apostolatos, Phys. Rev. D 89,124013 (2014), 1403.6243. 

[15] K. Yagi and N. Yunes, ArXiv e-prints (2015), 1503.02726. 

[16] J. M. Lattimer and M. Prakash, Apj 550, 426 (2001), astro-ph/0002232. 

[17] M. Urbanec, J. C. Miller, and Z. Stuchllk, Mon. Not. Roy. Astron. Soc. 433,1903 (2013), 1301.5925. 

[18] M. Baubock, E. Berti, D. Psaltis, and F. Ozel, Apj 777, 68 (2013), 1306.0569. 

[19] M. AlGendy and S. M. Morsink, ArXiv e-prints (2014), 1404.0609. 

[20] N. Andersson and K. D. Kokkotas, Mon. Not. Roy. Astron. Soc. 299,1059 (1998). 

[21] L. K. Tsui and P. T. Leung, Mon. Not. Roy. Astron. Soc. 357,1029 (2005), gr-qc/0412024. 

[22] E. Gaertig and K. D. Kokkotas, Phys. Rev. D 83, 064031 (2011). 

[23] D. D. Doneva, E. Gaertig, K. D. Kokkotas, and C. Kruger, Phys. Rev. D 88, 044052 (2013), 1305.7197. 

[24] T. Delsate, ArXiv e-prints (2015), 1504.07335. 

[25] P. Pani, ArXiv e-prints (2015), 1506.06050. 

[26] G. Pappas, ArXiv e-prints (2015), 1506.07225. 

[27] Y.-H. Sham, L.-M. Lin, and P. T. Leung, Apj 781, 66 (2014), 1312.1011. 

[28] B. Kleihaus, J. Kunz, and S. Mojica, ArXiv e-prints (2014), 1407.6884. 

[29] P. Pani and E. Berti, Phys. Rev. D 90, 024025 (2014), 1405.4547. 

[30] D. D. Doneva, S. S. Yazadjiev, K. V. Staykov, and K. D. Kokkotas, Phys. Rev. D 90,104021 (2014), 1408.1641. 

[31] G. Pappas and T. P. Sotiriou, Phys. Rev. D 91, 044011 (2015), 1412.3494. 

[32] G. Pappas and T. P. Sotiriou, ArXiv e-prints (2015), 1505.02882. 

[33] D. D. Doneva, S. S. Yazadjiev, N. Stergioulas, and K. D. Kokkotas, Phys. Rev. D 88, 084060 (2013). 

[34] N. Stergioulas, A. Bauswein, K. Zagkouris, and H.-T. Janka, Mon. Not. Roy. Astron. Soc. 418, 427 (2011). 

[35] G. Pappas and T. A. Apostolatos, Physical Review Letters 108, 231104 (2012). 

[36] J. Lattimer, Annu. Rev. Nucl. Part. Sci. 62, 485 (2012). 

[37] A. W. Steiner, J. M. Lattimer, and E. F. Brown, Apj 722, 33 (2010). 

[38] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts, and J. W. T. Hessels, Nature 467,1081 (2010). 

[39] J. Antoniadis, P. C. Freire, N. Wex, T. M. Tauris, R. S. Lynch, et al., Science 340, 6131 (2013). 



8 


[40] S. Yazadjiev, in Conference: (Non-)Universal Properties of Neutron Stars, ZAMS, Bremen (2015). 

[41] E. Barausse, C. Palenzuela, M. Ponce, and L. Lehner, Phys. Rev. D 87, 081506 (2013). 

[42] C. Palenzuela, E. Barausse, M. Ponce, and L. Lehner, Phys. Rev. D 89, 044024 (2014). 

[43] M. Shibata, K. Taniguchi, H. Okawa, and A. Buonanno, Phys. Rev. D 89, 084005 (2014). 

[44] K. Taniguchi, M. Shibata, and A. Buonanno, Phys. Rev. D 91, 024033 (2015), 1410.0738. 



